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Abstract 

In this paper we consider a heavy gluino to be the lightest super- 
symmetric particle [LSP]. We investigate the limits on the mass of 
a heavy gluino LSP, using the searches for excess events in the jets 
plus missing momentum channel in Run I. The neutral and charged 
R-hadrons, containing a heavy gluino LSP, have distinct signatures at 
the Fermilab Tevatron. The range of excluded gluino masses depends 
on whether the R-hadron is charged or neutral and the amount of 
energy deposited in the hadronic calorimeter. The latter depends on 
the energy loss per collision in the calorimeter and the number of col- 
lisions; where both quantities require a model for R-hadron- Nucleon 
scattering. We show how the excluded range of gluino mass depends 
on these parameters. We find that gluinos with mass in the range 
between ~ 35 GeV and ~ 115 GeV are excluded by CDF Run I data. 
Combined with previous results of Baer et al., which use LEP data 
to exclude the range 3 - 22~25 GeV, our result demonstrates that 



an allowed window for a heavy gluino with mass between 25 and 35 
GeV is quite robust. Finally we discuss the relevant differences of our 
analysis of Tevatron data to that of Baer et al. 
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1 INTRODUCTION 



Supersymmetric theories are the most promising candidates for new physics 
beyond the standard model. Since supersymmetry [SUSY] is not observed at 
low energies, it must be broken. Two interesting classes of models for SUSY 
breaking are minimal supergravity |l| (also known as the constrained minimal 
SUSY standard model [CMSSM] 0) and minimal gauge-mediated SUSY 
breaking [GMSB] |J. In both cases gauginos are assumed to have a common 
mass at the grand unification [GUT] scale ||. Since gaugino masses, at one 
loop, are proportional to their corresponding coupling constants, we naturally 
expect gluinos to be the heaviest gauginos at low energies. However, ref. 
[^, |6| has shown it is possible (and even natural) to build phenomenologically 
acceptable GMSB models with a heavy gluino as the LSP. [] In these GUT 
models, gaugino masses do NOT unify at the GUT scale. 

In GMSB models, supersymmetry is broken in a hidden sector and trans- 
mitted to the visible sector via messengers carrying standard model gauge 
interactions. The masses of gauginos arise via one loop diagrams containing 
messenger fields with appropriate quantum numbers. The minimal messenger 
sector includes color triplets, giving mass to gluinos and squarks, and weak 
doublets, giving mass to charginos and neutralinos as well as squarks and 
sleptons. These minimal messengers form complete five dimensional repre- 
sentations of SU(5); thus preserving GUT predictions. If the color triplet and 
weak doublet messengers are degenerate, then gaugino masses rhi,i = 1,2,3 
satisfy fhi/oii = constant at the messenger scale. If, on the other hand, the 
color triplet messengers are heavier than their weak doublet partners, gluino 
masses are suppressed. In the model of ref. || ||, it was argued that the 
latter possibility is "natural" in an SO (10) SUSY GUT with minimal mes- 
senger sector. In this case, the Higgs and messenger fields both belong to 
the 5 + 5 of SU(5) (or a 10 of SO (10)). Since they have identical quantum 
numbers it is natural for them to mix. Moreover the natural mass scale for 
real SU(5) representations is of order the GUT scale. Recall, in any SUSY 
GUT, Higgs triplet and doublet masses must be split in order to avoid rapid 
proton decay. It was shown in a simple example that the combination of 
Higgs-messenger mixing and doublet-triplet splitting leads to color triplet 

1 This is not the only model which may lead to a gluino LSP; see for example, the so- 
called O-II string model discussed in ref. [fij or another GMSB model discussed in ref. || . 
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Higgs and messengers with mass at the GUT scale (Mq) and lighter doublet 
messengers with mass at a scale M < Mq- The Higgs doublets remain mass- 
less. M is identified as the messenger scale. In this scenario the gluino is the 
LSP. Assuming a conserved R-parity, such a gluino is stable. 

A natural framework for the above scenario is obtained in SO (10) with 
the Dimopoulos-Wilczek mechanism |J for doublet-triplet splitting. This 
theory requires two 10s, the Higgs 10 dimensional multiplet and an auxiliary 
10 field. The Higgs and auxiliary fields mix; and if the auxiliary 10 feels SUSY 
breaking at tree level we have GMSB with the auxiliary 10 identified as the 
messengers 0]. Since M is both the messenger scale and the R symmetry 
breaking scale, the gluino mass is suppressed by (M/Mg) 2 compared to other 
gaugino masses. Note, in this theory, as shown in ref. ]TU|], M/Mq < 0.1 is 
also sufficient to suppress the baryon number violating nucleon decay rates. 

Given the possibility of having a gluino LSP it is important to consider 
the existing limits. [] In a much studied class of models |13| gauginos are 
massless at tree level obtaining their mass radiatively at low energies. In 
these radiative gaugino models the lightest neutralino is the LSP with a 
gluino NLSP (next to the lightest superparticle) with mass of order a couple 



of GeV or less. These models provide candidates for dark matter [fj] and 



UHECRons |L5|| , the particles responsible for the Ultra High Energy Cosmic 



Rays observed by cosmic ray shower detectors. There are two classes of 
experimental constraints on these models. The first class looks for the decay 
of the gluino into a quark- anti-quark pair and the LSP. Such a search by 



KTeV 1 16] has ruled out these models in the entire region of parameter space 



consistent with the dark matter and UHECRon solutions. The second class 
of constraints is valid even if the gluino were the LSP and thus stable. For 
example, beam dump experiments at Fermilab, as well as, searches for the 
decay T — > rjg + 7 rule out gluinos with mass in the range ~ 2 — 4 GeV. In 
addition, gluinos in the entire light gluino window 1 — 5 GeV are ruled out by 
searches for dijet events at Fermilab for squark masses roughly in the range 
150 — 600 GeV |17]]. Finally the analysis of the running of a s from m T to Mz 
combined with multijet angular correlations on the Z has lead the authors of 



IS] to claim that the entire light gluino window is ruled out. The analysis 



2 There are also cosmological constraints on a gluino LSP|Tl). These constraints, if 
correct, would require the gluino to be the NLSP and to decay, for example, into a gravitino 
(such as discussed in refs. || [l2|]) or into a neutralino as in 
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of the multijet angular variables however has large uncertainties as noted in 



IS]. If we neglect this part of the analysis, only using the running of a & 



then the light gluino window is ruled out only at 80% CL []18||. In summary, 



the light gluino window may still be viable for stable gluinos. 

In this paper, however, we only consider models with a heavy gluino LSP 
with mass greater than 5 GeV. In fact, as has been discussed in ||, all 
previous searches ruling out heavy gluinos in the MSSM (with mass greater 
than 5 GeV) use the jets plus missing momentum signature coming from 
the sequential decay of a gluino into q + q + the LSP (typically the lightest 
neutralino). However, if the gluino is the LSP and thus stable (assuming 
a conserved R-parity) then all previous limits on a heavy gluino, and most 
searches for SUSY, must be re-evaluated. Hence it is important to study a 
heavy gluino LSP (in the heavy gluino window). 

Current experimental data can greatly constrain these models [EL H], 
[ f20|j - |25|. An important new study of the heavy gluino LSP has been 



carried out by Baer, Cheung and Gunion [BCG] BCG have significantly 
constrained the heavy gluino window, using both OPAL and CDF data; their 
results are presented below. An important parameter in their analysis is P 
or 1 — P, the probability for a gluino to fragment into a charged or a neutral 
R-hadron.Q The same parameter, in their analysis, is used for the probability 
of the final state R-hadron in an R-hadron - Nucleon collision to be charged 
or neutral. 



Using OPAL data j^J for e+e~ ->• Z -»• N 2 N t with N 2 -> A x + q + q 
where the N 2 (N 1 ) neutralino is the NLSP (LSP) in the MSSM: This 
process was analyzed by OPAL for 2, 3 or 4 jets + missing momen- 
tum. Analyzed in terms of the process e + e~ — > qqgg BCG find that 
gluino masses in the range 3 - 22-25 GeV are ruled out. This result is 
independent of the probability P. 

Using CDF data for jets + missing momentum ffi'J: BCG find that 
gluinos with mass in the range 25 to 130 - 150 GeV are ruled out for all 
values of P in the range < P < 3/4. Combined with the OPAL data, 
this excludes a heavy gluino LSP in the range from 3 to 130 - 150 GeV. 
However for values of P > 3/4 and large hadronic scattering length, 
they find an allowed heavy gluino window from 23 - 50 GeV. 



3 An R-hadron denotes the color singlet bound state of a gluino with gluon or light 
quark constituents. 
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It is the goal of this paper to better understand the sensitivity in the BCG 
analysis to the probability parameter P and hence the existence of the heavy 
gluino window. Thus we re-analyse the CDF data using an entirely different 
model for R-hadron - Nucleon scattering. As in the analysis of BCG, we first 
assume the gluino LSP is heavy and that all other supersymmetric particles 
are substantially heavier than the gluino. We then check the sensitivity of 
our results to lighter squarks. Our results are found in sec. |j. The basic 
conclusion is that there exists a heavy gluino window in the mass range from 
25 - 35 GeV. 

The paper is organized as follows. We review the processes responsible 
for R-hadron production and fragmentation at the Tevatron in sec. [2.1| and 
for detection at CDF in sec. [272|. In particular, the Regge model used for 



R-hadron - Nucleon scattering is found in sec. |2.2.2| . In sec. |]we present the 
details of our analysis, i.e. evaluating missing momentum; ionization energy 
loss, and muon identification. Our results are presented in sec. £| and the 
conclusions are given in sec. [5]. 



2 A HEAVY GLUINO LSP AT THE TEVA- 
TRON 

2.1 Production and Fragmentation 

If gluinos are light enough, they can be produced in proton-antiproton col- 
lisions at the Tevatron. The relevant Feynman diagrams are shown in Fig. 

B 

Gluinos are color octets; hence gluinos cannot exist in isolation. They 
form colorless bound states such as R = gg (gluino-gluon bound state, an 
isospin singlet) or p = g(qq)(i=i) (gluino-quark-antiquark bound state with 
isospin 1) with (p + = (gud), p° = (guu — gdd) /y/2, p~ = (gdu)). In our 
analysis, we assume that these four states are the lightest bound states of 
the gluino. Together, they are called R-hadrons. As we see, two of these four 
R-hadrons are electrically charged and two are neutral. We also assume that 
any other bound state of the gluino, when produced, rapidly decays into one 
of these four R-hadrons. 

If nip > mn + m n then the decay p — > Rqtt will occur on a strong 
interaction time scale of order 10 -23 s. Otherwise p ± — > Rq^u or p° — > Roj 
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Figure 1: Diagrams representing gluino pair production in p — p collisions. 
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occur via weak or electromagnetic interactions. 
For p + we obtain the decay rate 



r(/?+ - J2VV) = f^xr^eV) (i) 



where the ratio 



r(p-e+e-) 15 a 2 M 2 R 



(2) 



AM = rrip± — m Ro , and M R and m M is the R-hadron and up quark mass, 
respectively. Using the experimental value 

T(p -> e + e~) = 0.00675 MeV (3) 

we find the lifetime t~ p + = 10~ 7 s for M R = 50 GeV and AM = 140 MeV or 
CTp± ~ 30 m. 

Similarly for p° we have 



lV^° + 7 )= ng^jl x r(p - e + e -) (4) 



and we obtain [] 



_ r(p°^7 + ^ + ^) _ 8 ^ (7r 2 -9)«2 AM 



r = 



x y —^-T7T-I (5) 



T(p -> e+e~) 45 tt a M|m2 



Using the same numerical values as above, we find the lifetime Tpo ~ 10 _14 s 
or cTpO ~ 0.0003 cm. Note for smaller values of AM the lifetime quickly 
grows. For AM = 5 MeV, we have crpa ~ 200 cm. 

Gluinos produced in a proton-antiproton collision fragment into one of 
these bound states. The probabilities of fragmenting into R , p + , p° and p~ 
are denoted by P R o, P p +, P p o and P p -, respectively with P R o + P p + + P p o + 
P p - = 1. If a particle lifetime is too short for the particle to multiple scatter 
in the detector, then the effective fragmentation probability for that particle 



is zero (see discussion of At in sec. j2.2.1| ). For example, if mp > m Ro + 
then we have P Ro = 1. On the other hand for mp < m Ro + m n and taking 
AM = 140 MeV and M = 50 GeV in the equation above we find ergo ~ 



4 For this calculation we use the analysis leading to equation (4) |2£| on the radiative 
decay J/* ^ jgg. 
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3 x 10 -4 cm giving P p o « 0. If all non-vanishing probabilities are then taken 
approximately equal we find the probability P (as in BCG) for producing a 
charged R-hadron given by P ~ 2/3.Q If, however, we take AM = 5 MeV 
instead, then all four probabilities are effectively non-zero and P ~ 1/2. In 
general, the effective fragmentation probabilities depend on AM, M and the 
relativistic 7 factor. In our analysis, however, we take these fragmentation 
probabilities to be free parameters and we test the sensitivity of our results 
to changes in them. 



2.2 CONSTRAINTS FROM CDF 

To constrain the mass of the heavy gluino LSP we use, following the analysis 
of Baer et al. [12], CDF data for jets+/> T ||27|| . This data has been used 



by CDF to put limits on gluino and squark masses. In the MSSM a heavy 
squark or gluino typically decays sequentially producing several jets plus 
missing momentum carried away by the LSP. In addition, in order to reduce 
the standard model background for the jets + fe T signature, CDF cuts all jets 
containing a charged lepton. 

In our case the lightest R-hadron is the LSP and since it is a hadron (in 
contrast to the CMSSM) it will deposit some of its energy in the hadronic 
calorimeter. In addition, the R-hadron's charge may fluctuate due to hadronic 
collisions. Thus as emphasized by BCG it is important to analyze in detail 
how an R-hadron is observed in the detector. In fact it was shown that some 
of the time a charged R-hadron may be identified as a muon. 

The important components of the CDF detector (for our analysis) are the 
central tracker, hadronic calorimeter, the inner and outer muon chambers 
and the un-instrumented piece of iron, located between the muon chambers. 
Consider how an R-hadron is observed. A gluino fragments into one of the 
four stable R-hadrons.f] If the R-hadron is charged, it leaves a track in the 
central tracker. Then the R-hadron enters the hadronic calorimeter (layers 
of iron) where it scatters off the nucleons a few times before it is stopped or, 
if not, it enters the inner muon chamber .[] After each hadronic collision, the 

5 Note, in our analysis these probabilities are only used for the initial fragmentation. 

6 Note, typically the R-hadron carries most of the gluino momentum ( ~ 99.4% for a 50 
GeV gluino); leaving little or none for associated pions or other hadrons. Thus a gluino 
jet often contains only the single R-hadron and nothing else. 

7 As discussed in ref. Q , it is reasonable to take the individual nucleon as the target 
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R-hadron can change into any other possible "stable" gluino bound state, i.e. 
any of the four R-hadrons with the probability determined by the appropriate 
scattering cross section. A charged R-hadron also deposits ionization energy 
in the hadronic calorimeter. Note, if the R-hadron exiting the hadronic 
calorimeter is charged, it is detected in the muon chambers. In addition, 
between the first and second chambers it can lose energy and also the charge 
of the R-hadron can change due to an hadronic interaction inside the un- 
instrumented iron. The energy loss will again be both in the form of hadronic 
energy loss and ionization (if it is charged). However, this energy loss is not 
detected. 

In order to determine the amount of energy deposited in the hadronic 
calorimeter and to model the charge fluctuations in the detector, we need 
a model for R-hadron scattering. The magnitude of the total cross section 
ot determines the scattering length \ T oc 1/<Jt- The angular dependence of 
the differential scattering cross section determines the mean energy loss per 
collision AE. Finally the differential cross sections determine the probability 
for charge fluctuations. In these details reside the main differences between 
our analysis and that of BCG. 



2.2.1 The total cross section and the interaction length At- 

A priori, we do not know the magnitude of <jt- However we estimate it using 
the two gluon exchange model for the total cross section of R-hadron-nucleon 
collisions developed in ref . |JD| and used in [JT2] . This analysis suggests that 



the ratio 

X T (R) _ a T {irN) _ C F < r 2 w > 

A T (vr) a T (RN) K C A )K <r 2 R > } 1 ' 

where Cp = 4/3(C*a = 3) refers to the quadratic casimir of SU(3) co ; or in 
the fundamental (adjoint) representations and < r 2 > is the transverse size- 
squared of the particle. 

The first factor is due to the octet nature of R constituents. The reduced 
constituent mass of the particle determines < r 2 >. In the case of a pion 
we have < r 2 >oc 4/m 2 and in the case of R° we get < r 2 R0 >oc l/m 2 
(for nig ^> trig) , where m q and m g are the light quark and gluon masses, 

for values of 1 > \t\ > 0.01 GeV 2 . For smaller values of |t| the entire nucleus would be 
treated coherently and for larger \t\ individual quarks would be the scatterers. 
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respectively. BCG assume equal values for m q and m g and estimate 



'_R°Af (9/16)0"^ which translates into Xt(R°) ~ (16/9)At(vt). Noting that 
At(vt) ~ 11 cm, they conclude that Xt(R°) ~ 19 cm is a reasonable estimate. 
However the constituent quark mass is roughly speaking (l/3)m proton ~ 330 
MeV and the constituent gluon mass is (l/2)m g i ue b a u ~ 750 MeV. Hence 
it might be more realistic to use the relation m g ~ 2m q which changes the 
result to Xt(R°) ~ 78 cm. In this analysis we take At = 19 cm and Ay = 38 
cm, which are the values used by BCG 

It will be seen later that larger values of Xt{R°) can open an allowed 
region in the gluino mass window. Although we have adopted the estimates 
of BCG, we believe that choosing a larger value for At is NOT an extreme 
limit, but is in fact preferred. We also assume that all R-hadrons have the 
same At- We rescale all the differential cross sections to get the desired Al- 
and keep their ratios unchanged. Note, the ratios are determined by the 
Regge analysis discussed in the following section. 

The average number of hadronic collisions inside the hadronic calorimeter 
depends on the magnitude of At- For At = 19 cm, the R-hadrons on average 
undergo 6 hadronic collisions in the hadronic calorimeter and 2 hadronic 
collisions in the un-instrumented iron between the muon chambers. 



2.2.2 Regge scattering model: AE and charge fluctuation 

The differential cross sections for the four R-hadrons determine both AE, the 
energy loss per hadronic collision and the probability for charge fluctuation. 
In this analysis we use single particle inclusive Regge cross sections with 
Pomeron and p Regge exchanges. Recall that Regge exchange is by far the 
dominant contribution to high energy hadronic scattering cross sections. 

Let us briefly review Regge theory. The Regge pole model connects the 
two classes of phenomena: 

• classification of hadrons, and 

• high energy scattering of hadrons. 

Hadronic bound states or resonances with identical quantum numbers, 
except spin and mass, are correlated by Regge trajectories where they ap- 
pear with spin J differing by two units. There is another relevant quan- 
tum number in Regge theory known as the signature r. For mesons, signa- 
ture is defined to be r = (— 1) J with the corresponding particle sequences, 
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J = 0, 2, 4, (for r = +1) and J = 1,3, 5, (for r = —1). As proposed 



by Chew and Frautschi f31fl , particles are classified on Regge trajectories in 
plots of their spin J versus the square of their mass M 2 . The trajectory a, 
interpolates between the particles such that Re[a(Mj)] = J.f\ The trajec- 
tory a(t) is assumed to be an analytic function of t and experimentally they 
are, to a good approximation, straight lines. The two body scatterings are 
dominated by the exchange of one or more of these trajectories. To see how 
a Regge pole exchange amplitude could arise in a field theory, consider the 
Van Hove-Durand model for the scattering of A + B ^ A + B with A and 
B spinless particles. A and B interact via exchange of a meson with spin J 
in the t-channel. The amplitude is given by 

Ms t ) - 9 2 j[-s/s,] J 

in the large s limit. If the exchanged particle is a member of a Regge tra- 
jectory having an infinite series of mesons with spins J = 1,3, 5,..., the scat- 
tering amplitude can be expressed as a summation of the single exchanges of 
all the mesons on the trajectory. 

where sq is an energy scale factor. Assuming a universal coupling and a linear 
Regge trajectory Mj = /i 2 (J — a), the series sum can be expressed in closed 
form. The result in the limit of large s is 

2fi z sm(7ra) s 

which is the form for the exchange of Regge poles with odd signature with 
a trajectory a(t) = a + tf /i 2 satisfying a(Mj) = J. The trajectory a(t), for 
t < describing scattering in the s channel is a smooth continuation of the 
Regge trajectory for t > describing scattering due to resonance exchange 
exchange in the t channel. In general, the form of a Regge amplitude can be 
more complicated if the incoming particles carry spin or charge. 

8 We only talk about meson trajectories here; baryon trajectories are similarly discussed 
inref. EJl. 
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Figure 2: Diagrams showing an inclusive process a + fe— >c + X in terms of 
a Regge pole a(t) in the t-channel 

2.2.3 INCLUSIVE SCATTERING, TRIPLE REGGEON FORM 

The single particle inclusive scattering process a + b —> c + X can also be 
expressed in terms of Regge exchanges, a, h and c are definite particles and 
X represents everything else with total invariant mass Mx- In the limit 

M x large, s/M x > 1 and s > t (7) 

the near forward two-body inclusive scattering process can be described by 
the t-channel exchange of a Reggeon. Although single particle inclusive scat- 
tering represents only ~ 20% of the total hadronic cross section we will use 
this model for R-hadron scattering in our analysis. It determines the energy 
loss AE as well as the probability for charge fluctuation. This process is 
represented in Fig. |2] . 

The amplitude for the process in Fig. | is given by 



A{s,t) <xp b£ (t)£(t)s a V. 



(8) 



Using a generalized optical theorem due to Mueller JHJ as shown in Fig. 
|3| the inclusive cross section a + b c + X is related to the discontinuity in 
the forward three-body amplitude a + b + c-^a + b + c. In the Regge limit 
(eqn. [?]) this gives the triple Regge cross section (see Fig. |j) 



d 2 a 
dxdt 



i,j,k 



(*)( 



Mi' 



\ai(t)+aj(t) 



(m; 



|£»*(0) 



(9) 
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a. 



a 



Figure 3: Using a generalized optical theorem due to Mueller [^] we relate 
the single particle inclusive cross section a + b^c + X to the discontinuity 
in the forward three-body amplitude a + b + c^a + b + c 



with x = 1 — ^ JL . 

s 

For simplicity we consider two Regge trajectories, an isosinglet Pomeron 
trajectory P and an isovector Reggeon trajectory p. We then write all the 
possible triple couplings of the Pomeron, Reggeon and R-hadrons as follows 

g RRP R l Q RlP + g-ppp ~p\p\P + g R p P Ropipi + g- P p P e m p\p)p k , ( 10) 

where the indices 1 and 2 distinguish between the outgoing and incoming 
particles in a vertex. 



Using 



we find 



P± 



PiPi = 

PiPi = 
ttjkp\p)pk 



Pi ± ip2 
V2 



P± 



Pi ± ip2 
V2 



P+P- + P~P+ + P3P3, 

P\p- + p-p 2 + + pIpI 



(ii) 



(12) 



K+P\p-Pz 
-p\~p\p- 



p-p 2 +P3 + P-pIp+ 



pXp\p- 



We thus obtain the following vertices and the relevant couplings given by 

iZo^o^P with g RRP , (13) 
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Figure 4: The triple Reggeon scattering diagram describing the diffractive 
process a + b — > c + X 

p+p-P, p-p+P and p 3 p 3 P with g~ p ~ pP , 

R p-p + , Rop+p- and R0P3P3 with g Rpp , 

+P+P-P3, -P-P+P3, +P-P3P+, 
-p+p 3 p_, -p 3 p + p_ and + p 3 p_p+ with %g ppp . 

which lead to the following diffractive scattering patterns, 

RqN — >RqX ) - Gppp, 

R N-^p~( + _ 3 )X, p (+ _ 3 )iV— >_R X =>• Gppp, 
P+N-+P+X, p^N^p-X G 2 PPP} G 2 ppP , 
p 3 N^p 3 X Gppp, 
p + N^p-X, p^N^p + X = 0, 
p + N-*p 3 X, p 3 N^p + X, p_N^p 3 X, p 3 N^p_X G 2 ppP . 
Using eq. ^| and the Pomeron and Reggeon trajectories 

op(i) = 1 + a P t, a p (t) = a + (3t, (14) 
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we obtain the desired differential scattering cross section formulas for the 
triple Pomeron contribution, 



dM^dt ~ 2GpPp{t){ W ( M 



2G PPP (t)( T ^Y^( JI -), (15) 



and for the Reggeon-Reggeon-Pomeron contribution, 



These results are used to perform the analysis in this paper.| The calculations 
are based on the assumption that 

Gp PP (t) = G PPP (t) = Gp PP (t), G pp p(t) = G ppP (t) = G 2 ppP (t). (17) 



Inspired by ref. [34 , we take 



G P pp{t)/Gppp{0) = 0.88 e 3 - 94 * + 0.12 e 1 ' 12 * (18) 
Gppp(t)/Gp pP (0) = 0.85 e 7 - 26 * + 0.15 e" 1 - 83 *. 

and 

a P = 0.36 (19) 
a = 0.5 
/3= 1.0 . 

We also vary the relative size of Gppp(t) and G pp p{t) to check the sensi- 
tivity of our results to these couplings. 

Before continuing, let us briefly discuss the scattering analysis of BCG 



1^]. These authors assume that when an R-hadron scatters on a nucleon 
the light brown muck (quarks and gluons bound to the gluino) are stripped 
off. The bare gluino then re-fragments into a charged R-hadron with the 
fragmentation probability P and into a neutral R-hadron with probability 
1 — P. Thus for BCG, the probability P plays two independent roles: 

9 For simplicity we have ignored the triple Reggeon contribution. We also use the Regge 
cross-sections for all physical values of s, t and Mx; in particular, extending outside the 
Regge limit (eqn. 2). 
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• fragmentation probability, and 

• rescattering probability. 

In our analysis, however, Regge scattering cross sections allow R-hadrons 
to change their identity with velocity dependent probabilities in each hadronic 
collision. We thus separate the independent phenomena of fragmentation and 
rescattering. 

In addition, to calculate the energy of an R-hadron after an R-hadron - 
nucleon collision, BCG use either a constant differential cross section 
which vanishes for t > 1 GeV or a triple Pomeron scattering formula. We on 
the other hand use the differential Regge cross sections for this as well. 



3 THE ANALYSIS 

Consider the CDF data for jets+ j6 T [27]. We compare our Monte Carlo 



simulation with this data and thus we use the same cuts as BCG, which were 
designed to duplicate the experimental procedures of ref. [p^j . The cuts are 
listed as follows. 

1. No (isolated) leptons with E T > 10 GeV. 

2. p T > 60 GeV. 

3. There are three or more jets with \r]j e t\ < 2 and Et > 15 GeV, using a 
coalescence cone size of AR = 0.5. 

4. Azimuthal separation requirements: A0(j6 T , ji) < 160° and 
A0(j6 T , j(E T > 20GeV)) > 30°. 

For our monte carlo analysis we generated events using (SPYTHIA, A 
Silver symmetric Extension of PYTHIA 5.7 |3"5|), which has been modified 
by S. Mrenna and K. Tobe to accommodate a gluino LSP. We have written 
a toy calorimeter code extending out to \r]j et \ < 4. The cell size of the 
calorimeter is set to At] x A<p = 0.1 x 0.1. The hadronic resolution is also 
taken to be 70%/\/^. The missing transverse momentum cut (2) and the 
jet-number cut (3) are the strongest cuts, jl is the most energetic jet in an 
event and j(Et > 20GeV) is any jet with a transverse energy larger than 20 
GeV. Cut (4) eliminates the QCD jet mis-measurement backgrounds. These 
backgrounds occur due to the uncertainty in the energy measurement of 
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very energetic jets leading to large missing momentum in either the same or 
opposite direction as their momenta. 



3.1 Identifying missing momentum 

From the total energy deposited in the hadronic calorimeter and the pseudo- 
rapidity T) of the corresponding jet, one can determine the measured trans- 
verse energy of the jet. We assume (like CDF) that the jets are massless. 
Hence the magnitude of the transverse momentum is the same as the trans- 
verse energy. 

If we add up the transverse momenta of all jets present in a single detected 
event, we should get zero due to the conservation of transverse momentum. In 
practice this does not happen. For one reason, not all the generated particles 
are energetic enough to be detected by a calorimeter cell as part of a jet (this 
effect is parameterized in pythia). The second and very important reason 
is that an R-hadron does not deposit all its energy inside the calorimeter; 
resulting in some missing momentum. For example, (1) only the kinetic 
energy of an R-hadron is visible in a detector and not its rest mass and (2) 
the energy loss per collision decreases as the R-hadron mass increases. There 
is also an enhancement in the ionization energy loss for R-hadrons considered 



later in Sec. p.2[ The total missing transverse momentum fa T is defined as 
minus the sum of the transverse momenta of all jets in each event. A large 
missing transverse momentum is a signature for a gluino LSP. 

Spythia gives us information which is directly converted into px for all 
standard model particles. For R-hadrons on the other hand the output of 
Spythia is input for our hadronic calorimeter code. It is in this code that At 
and Regge cross sections are used to determine the energy deposited in the 
hadronic calorimeter. As discussed earlier, R-hadrons suffer both hadronic 
and ionization energy loss in the hadronic calorimeter. Hadronic energy losses 
are due to collisions of R-hadrons with nuclei. We assume an R-hadron - 
nucleon collision occurs in the hadronic calorimeter after each distance At 
in the detector. We use the single particle inclusive differential Regge cross 
sections, discussed in sec. |2.2.3| , to obtain the energy and charge of the 



R-hadrons after each hadronic collision. These are computed randomly on 
an event by event basis using a probability distribution which is weighted by 
the differential cross sections. Each R-hadron can scatter into any of the four 
R-hadrons after colliding with a nucleon. Therefore there are four separate 
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processes which can occur. To decide which one occurs, we compute the total 
cross section for each process. The ratio of the total cross sections determines 
the probability of each process. A random number is generated which decides 
which process occurs based on these probabilities. In the next step we use the 
corresponding differential cross section to determine the energy of the out- 
going R-hadron. The energy of an out-going R-hadron is a function of the 
momentum transfer t and the invariant mass Mx of the inclusive process. 
The probability density for an R-hadron with given t, Mx and center of 
momentum energy y/s is d f^ x \ t,M x - A (t, M x ) pair is therefore generated, 
weighted by d f^ \t,M x t° compute the energy of the out-going R-hadron. 



3.2 Ionization energy loss 

To calculate the ionization energy deposit, we use the standard dE/dx for- 
mula from ref . []36| . A 50 GeV muon beam, as measured in the central track- 
ing chambers by CDF, deposits 2 GeV of energy in the hadronic calorimeter. 
If we compute the energy loss of the muon beam using the standard dE/dx 
formula for a muon passing through an iron detector, we find it deposits only 



1.3 GeV. Ref. |T2| thus defines the ratio of the measured calorimeter ion- 
ization loss to the actual dE/dx loss from ionization, E ion i zation , as r = 1.6. 
Finally, the visible energy of a jet (that measured in the calorimeter) is the 
energy we assign to the jets when calculating fa T . It is given by the total 
energy loss E total = rE ionization + E hadronic . As mentioned earlier all jets are 
massless and we identify the magnitude of the 3-momentum of a jet with 
Etotai and choose its original direction as the direction of its momentum 3- 
vector. In this way, we can add up the 3-vector momenta of all jets present 
in a single event and determine the total fa T . 



3.3 R-hadrons identified as muons 

Some of the R-hadrons which pass through the CDF detector are identified 
as muons. An R-hadron jet is declared to be muon-like, if it is detected to 
be charged in the central tracker and in at least one of the muon chambers, 
covering only a portion of the available pseudo-rapidity region. In addition, 
it should have a momentum larger than 10 GeV in the central tracker and 
also should not deposit more than 6 GeV inside the hadronic calorimeter 
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as measured by the detector. Each event containing a muon-like jet is dis- 



carded p7| . In some cases a large number of events are discarded because 
they contain muon-like jets. 

It is thus important to understand under what circumstances an R-hadron 
jet will be muon-like. Neutral R-hadron jets will not be muon-like since the 
jets need to be charged both in the central tracker and in one of the muon 
chambers. On the other hand, if a jet remains charged for most of its passage 
through the hadronic calorimeter, it is likely that it will deposit too much 
energy (more than 6 GeV) and thus will not be a muon candidate. Hence an 
ideal muon-like jet is one that is charged in the central tracker and one of the 
muon chambers and is neutral for most of its passage through the hadronic 
calorimeter. From the above discussion, we understand that increasing P p + 
and P p - will increase the number of R-hadron jets that are identified as muon 
jets. However if appp ^> a ppP (i.e. Pomeron scattering dominates) then 
there is no chance for charge fluctuations and typically the charged R-hadron 
deposits too much energy in the hadronic calorimeter. Finally, if At increases 
then there are fewer hadronic collisions and thus less hadronic energy loss in 
the calorimeter. In addition, as a consequence of fewer hadronic collisions, 
the R-hadron retains its velocity. This results in less ionization energy loss 
since a fast moving charged particle loses less energy due to ionization than 
a slower particle, if its velocity is below minimum ionizing which is the case 
most of the time here. This tends to increase the number of muon-like events. 



4 THE RESULTS 

To determine the limits on the mass of the gluino LSP, we compare the cross 
section for gluino pair production, after applying the cuts, with the standard 
model background as reported in the analysis by CDF . They report 



a background rate of 28.7 events Q for an integrated luminosity of 19 pb~ l 
at the Fermilab center of mass energy of y/s = 1.8 TeV corresponding to 
op = 1.51 pb. A 1.96cr background fluctuation corresponds to an allowed 
signal of a s ~ 553 fb, which is also equivalent to a s ~ 614 fb with an 
integrated luminosity of 0.1 fb" 1 at 5a level. This is the bound used by 



10 This background rate is determined using a Monte Carlo calculation. The measured 
background is actually 36 events. 
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Table 1: Labels for the different parameter choices used in the figures with 
At (interaction length), the ratios of triple Regge vertices and the fragmen- 
tation probability vector P = (Pro, P p +, P p o, P p -). 



Labels 


Parameter values 


11, 12 


Xt = 19cm, At = 38cm 


gi, g2 


appp ~ cf pp p, oppp ^> a p P p 


pi, p2, p3 


P= (1/4, 1/4, 1/4, 1/4), P = (0,1, 0,0), P= (1,0, 0,0) 



BCG ||12|| . We also use the latter cross section and do our analysis for an 
integrated luminosity of 0.1 /if" 1 . 

If the gluino signal, after cuts, is larger than a s ~ 614 fb with an inte- 
grated luminosity of 0.1 fb' 1 then those specific values of the gluino mass, 
At, fragmentation probabilities and Reggeon couplings, are excluded. We 
first analyze this parameter space in the limit that all other superparticles 
are too heavy to be produced at the Tevatron. We then determine the sen- 
sitivity of our results to lighter squarks by lowering the squark masses. At, 
Opppjo-ppp and the fragmentation probability vector P = (Pro, P p +, P p o, P p ~) 
are also varied to check the sensitivity of our results to these parameters. 
The labels assigned to the different parameter choices discussed below are 
given in table |I[ 

Fig. H corresponds to At = 19 cm, appp ~ a pp p and Pro = P p + = P p o = 
P p - = 1/4. Using our labeling system, this is called a (llglpl) scenario. As 
seen in this case, we are able to exclude gluino masses from 30 GeV up to 
130 GeV at 95% confidence level. The large signal cross section at 30 GeV 
suggests that the gluino mass exclusion can be pushed to even lower values. 
Considering the fact that BCG have excluded gluinos with mass in the range 
3 to 22 — 25 GeV using LEP data, the gluino window is closed from 3 to 130 
GeV for these scattering parameters. It will become evident that the most 
important parameter is At- 

Fig. |6| corresponds to (llg2pl). In this case an R-hadron which scatters 
off the nucleon will almost always retain its identity. The charged R-hadrons, 
produced 50% of the time in this case, will always remain charged as they pass 
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through the calorimeter. In (llglpl), it is possible for a gluino to fragment 
into a charged R-hadron but become neutral in the muon chambers. This 
will not happen for the case of (llg2pl). We might expect an increase in 
the number of muon-like jets in (llg2pl) as compared to (llglpl). On the 
other hand, there is a large ionization energy deposit for the charged R- 
hadrons in (llg2pl) because they remain charged all the way through the 
calorimeter. This can suppress the number of muon-like jets due to the fact 
that it would be less likely for them to deposit less that 6 GeV of transverse 
energy in the calorimeter. The second effect dominates over the first and we 
see a decrease in the total number of muon-like jets in (llg2pl) compared to 
(llglpl). Looking at Fig. [| and Fig. ^ we can see this effect which shows 
itself as an increase in the signal cross section, mostly at large rrig, in Fig. ^. 
In (llg2pl) the gluino mass is excluded up to 140 GeV. 

Fig. is labeled (llglp2) where, in contrast to (llglpl), all gluinos ini- 
tially fragment into charged R-hadrons. This greatly increases the chance of 
detecting a muon-like jet and results in a smaller signal cross section par- 
ticularly at larger masses. At smaller gluino mass, the missing momentum 
signal is reduced. This is a consequence of two effects; for lighter gluinos (1) 
the hadronic energy loss per collision is larger and (2) the unmeasured mass 
energy is smaller. Hence only a small number of events survive to be ana- 
lyzed for the muon cut at small masses and the difference due to the muon 
jet cut is not significant between (llglp2) and (llglpl) at small gluino mass. 
Overall, we see a smaller signal cross section in Fig. [7] compared to (llglpl) 
in Fig. |!| The result is that we only exclude gluino masses up to ~ 125 GeV. 

(Ilglp3) is opposite to (llglp2) where now all gluinos initially fragment 
into neutral R-hadrons. Hence there are no muon-like jets present in this 
scenario. We expect larger cross sections as seen in Fig. || However as in 
(llglp2), we do not see a significant increase in the cross section at small 
gluino mass. We now exclude gluino masses up to ~ 135 GeV. 

In all the above cases, the signal cross section at small gluino mass is 
similar for the reason explained previously. We therefore find, for Ay = 19 
cm, gluino masses are excluded, when combining LEP (T2| and CDF data, 
from 3 to at least 120 GeV. 

We now consider the analysis for Ay = 38 cm. We consider the (I2glpl), 
(I2glp2) and (I2g2p3) scenarios. 

The case (I2glpl) is plotted in Fig. || Larger At implies fewer R-hadron 
collisions with nuclei. This means that they generally move faster inside 
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Figure 5: The cross sections after cuts in the jets+ p T channel is compared 
to the 5<T level for L = 0.1 fb^ 1 which is roughly the same as the 95% CL for 
L — 19 pb^ 1 at yfs = 1.8 TeV. = 19 cm and the Pomeron contribution 
is comparable with the Reggeon contribution to the total scattering cross 
section. The gluino has initially fragmented into the neutral and charged 
R-hadrons with equal probabilities. 
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Figure 6: The cross sections after cuts in the jets+ p T channel is compared 
to the 5<T level for L = 0.1 fb^ 1 which is roughly the same as the 95% CL for 
L — 19 pb^ 1 at yHi = 1.8 TeV. Xt = 19 cm and the Pomeron contribution 
is much greater than the Reggeon contribution to the total scattering cross 
section. The gluino has initially fragmented into the neutral and charged 
R-hadrons with equal probabilities. 
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Figure 7: The cross sections after cuts in the jets+ fa T channel is compared 
to the 5a level for L — 0.1 fb^ 1 which is roughly the same as the 95% CL for 
L — 19 pb^ 1 at yfs = 1.8 TeV. At = 19 cm and the Pomeron contribution 
is comparable with the Reggeon contribution to the total scattering cross 
section. The gluino has initially fragmented only into the charged R-hadrons. 
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Figure 8: The cross sections after cuts in the jets+ fe T channel is compared 
to the 5a level for L = 0.1 fb^ 1 which is roughly the same as the 95% CL for 
L — 19 pb^ 1 at yt~s = 1.8 TeV. At = 19 cm and the Pomeron contribution 
is comparable with the Reggeon contribution to the total scattering cross 
section. The gluino has initially fragmented only into the neutral R-hadrons. 
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Figure 9: The cross sections after cuts in the jets+ p T channel is compared 
to the 5<T level for L = 0.1 fb^ 1 which is roughly the same as the 95% CL for 
L — 19 pb^ 1 at yHi =1.8 TeV. Xt = 38 cm and the Pomeron contribution 
is comparable with the Reggeon contribution to the total scattering cross 
section. The gluino has initially fragmented into the neutral and charged 
R-hadrons with equal probabilities. 
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Figure 10: The cross sections after cuts in the jets+ fa T channel is compared 
to the 5a level for L — 0.1 fb^ 1 which is roughly the same as the 95% CL for 
L — 19 pb^ 1 at yfs = 1.8 TeV. At = 38 cm and the Pomeron contribution 
is comparable with the Reggeon contribution to the total scattering cross 
section. The gluino has initially fragmented only into the charged R-hadrons. 
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Figure 11: The cross sections after cuts in the jets+ fa T channel is compared 
to the 5a level for L = 0.1 fb^ 1 which is roughly the same as the 95% CL for 
L = 19 pb^ 1 at ^/s = 1.8 TeV. At = 38 cm and the Pomeron contribution 
is much greater than the Reggeon contribution to the total scattering cross 
section. The gluino has initially fragmented only into the neutral R-hadrons. 
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the hadronic calorimeter, because they are slowed down less by hadronic 
collisions with the nuclei. The fast moving charged particles lose less energy 
due to ionization than slower particles. Thus there is less energy deposited in 
the (I2glpl) scenario than (llglpl). This significantly increases the number 
of muon-like jets as compared to (llglpl) because the jets are more likely 
to deposit less than 6 GeV transverse energy and they are identified muon- 
like if they satisfy the other requirements. The gluino mass is excluded up 
to slightly lower than 120 GeV. The signal rate also decreases substantially 
below 50 GeV. We observed only one event at 35 GeV. We thus conclude 
that we can not exclude a gluino mass below 35 GeV in our analysis. We 
therefore suggest that there is an open window between ~ 22 — 25 GeV and 
35 GeV requiring further exploration. 

In the (I2glp2) scenario (see Fig. [10]) we observe smaller signal cross 
sections than (I2glpl) for the same reason that (llglp2) has smaller signal 
cross sections than (llglpl). This is due to the greater number of events 
that contain muon-like jets. The gluino mass is excluded up to ~ 115 GeV. 
At small gluino mass we should not expect a significant difference between 
(12glpl) and (I2glp2) as mentioned before. Therefore, there is again an open 
window between ~ 22 — 25 GeV and 35 GeV in this case. 

The (I2g2p3) scenario is plotted in Fig. |TT|. There are no muon-like events 
in this scenario and we expect a larger signal rate compared to (I2glpl). We 
find events even at a gluino mass of 35 GeV and thus there is no open window 
below 35 GeV. The Gluino mass is excluded up to ~ 125 GeV. 



4.1 Lowering the squark masses 

To examine the effect of squarks on the results, we decrease the squark 
masses. A colliding quark and anti-quark make a virtual gluon in the s- 
channel which turns into a gluino pair. They can also exchange a squark 
in the t-channel and produce a pair of gluinos. When all squark masses are 
large, the former process dominates. Upon lowering the squark masses, the 
negative interference of the latter process decreases the gluino pair produc- 
tion cross section. However, the gluino production cross section due to the 
process gg — > gg remains unchanged. When squark masses are sufficiently 
decreased then squarks can be produced directly in proton-antiproton colli- 
sions. We therefore have events containing a gluino and a squark jet. The 
gluino-squark production cross section, however, is smaller than the gluino 
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pair production cross section by a factor of approximately 100 at large gluino 
mass (around 140 GeV) and a factor of 1000 at smaller gluino mass (around 
60 GeV). At large gluino masses, the signal rate increases substantially due 
to the squark production which enhances the jet signal rate and more than 
compensates the decrease in the total production cross section of gluino pairs. 
This however is not the case at smaller gluino masses. In fact the signal rate 
decreases upon lowering the squark masses (we have considered squarks as 
light as 450 GeV). At small gluino mass, squarks are produced at a relatively 
smaller rate than in the case of large gluino mass and the enhancement due 
to the squark signals are much less. The gluino pair production cross section 
decreases as mentioned before and the overall effect is that a lower signal rate 
is observed. Recalling that we found an allowed window for gluinos between 
~ 22 — 25 GeV and 35 GeV in the case of = 38 cm and equal probabilities 
for producing the four R-hadrons, the decrease in the signal rate has the ef- 
fect of enlarging the heavy gluino window, pushing it to slightly higher values 
of the gluino mass. 



5 CONCLUSIONS 

In this paper we considered a heavy gluino to be the lightest supersymmetric 
particle [LSP]. We used CDF Run I data, for excess events in the jets plus 
missing momentum channel, to place limits on the mass of a heavy gluino 
LSP. We find that gluinos with mass in the range between ~ 35 GeV and 
~ 115 GeV are excluded. Combined with the previous results of Baer et 



al. [BCG] 0, which use LEP data to exclude the range 3 - 22~25 GeV, 
our result demonstrates that an allowed window for a heavy gluino with mass 
between 25 and 35 GeV is quite robust. 

We compared our analysis to that of BCG. They also found an open win- 
dow, however, for seemingly unphysical values of the fragmentation proba- 
bility parameter P > 3/4. In their analysis, the bare gluino fragments into 
a charged R-hadron with the fragmentation probability P and into a neutral 
R-hadron with probability 1 — P. These authors also assume that when an 
R-hadron scatters on a nucleon the light brown muck (quarks and gluons 
bound to the gluino) are stripped off. The gluino then re-fragments with the 
same probability P. Thus for BCG, the probability P plays two independent 
roles: 
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• fragmentation probability, and 

• rescattering probability. 

In our analysis, however, Regge scattering cross sections allow R-hadrons to 
change their identity with velocity dependent probabilities in each hadronic 
collision. We thus separate the independent phenomena of fragmentation 
and rescattering. As a result, we find a heavy gluino window with physical 
fragmentation probabilities discussed in sec. ^TT| . 

Both our results and those of BCG are sensitive to the R-hadron - Nucleon 
scattering length At- The heavy gluino window only exists for At = 38 cm. 
For the smaller value of At = 19 cm, the window is closed; consistent with the 
results of BCG. Although we do not know a priori the value of the scattering 
length, we argued that larger values are preferred. 

Note, as At is increased, at some point the R-hadron will no longer have 
hadronic collisions in the calorimeter. The neutral R-hadron will escape the 
detector, while will only undergo ionization energy loss; behaving just 
like a heavy charged lepton. We do not expect the heavy gluino window 
between 25 and 35 GeV (discussed here) to be significantly affected. On the 
other hand, one might be concerned that searches for charged stable massive 
particles would now provide significant contraints on a heavy gluino in this 
limit. However BCG have performed an analysis of charged stable massive 
particles (in the context of the heavy gluino LSP) using CDF data [J37|J . 
They find that no limit can be set for a gluino mass less than 50 GeV. 

In addition, an open window requires significant R-hadron charge fluctu- 
ations. We showed that in the case of purely neutral R-hadrons, gluinos are 
excluded with mass between 3 and 120 GeV (see Fig. |PTD . 

Charge fluctuations are relevant for the analysis, since the CDF data cuts 
events with isolated leptons, in particular muons with > 10 GeV. A track 
is identified as a muon if it is charged in the central tracker and in at least one 
of the muon chambers. It's energy must also be more than 10 GeV, as mea- 
sured in the central tracker, but less than 6 GeV in the hadronic calorimeter. 
Thus R-hadrons can be identified as muons if they initially fragment charged 
and then end up in at least one of the muon chambers charged. However in 
order to deposit less than 6 GeV in the hadronic calorimeter, it is important 
to have a longer scattering length (fewer hadronic collisions) and significant 
charge fluctuations (less ionization energy loss). 

It is clear from our analysis that gluinos with mass in the range 25 - 
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35 GeV (and At = 38 cm) are cut from the CDF sample by the muon cut 
and even more so by the 60 GeV missing momentum cuts. This greatly 
diminished CDF's sensitivity to the heavy gluino LSP in this mass range. It 
may be possible to re-analyze the Run I data with revised cuts which enhance 
the signal to background ratio. Such an analysis is however beyond the scope 
of this paper. 

Consider alternate methods for finding the heavy gluino. We note that 
of order 10 6 R-hadrons, with mass in the allowed window, were produced 
at CDF in Run I. Assume only a few percent are stopped in the hadronic 
calorimeter. These will be absorbed into the iron nuclei and thus a small 
sample of the calorimeter can be tested in a mass spectrometer to search for 
heavy isotopes of iron. It is also possible for the gluino to be the next-to- 
lightest superparticle and decay into a gravitino and a gluon with lifetimes on 
the order of 100 years M . This possibility would not affect any of the analysis 
carried out in this paper{£] On the other hand, this possibility would result 
in a spectacular heavy gluino signal. The stopped R-hadrons would decay 
into a gravitino plus hadrons with visible hadronic energy niR — mGravitina 
which could be of order 10 GeV. Clearly one can now look for events in the 
detector, when the accelerator is OFF. 

In conclusion, it is important to search for a heavy gluino LSP with 
mass in the allowed range 25 - 35 GeV. After all, this may be where the 
elusive SUSY is hidden. In addition, as discussed by Albuquerque et al. 
these particles are the prime particle physics candidates for the so-called 



UHECRons [15 , responsible for the observed ultra high energy cosmic ray 
showers. 
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